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Abstract 

Kinematic edges in cascade decays provide a probe of the masses of new particles. In some new 
physics scenarios the decay chain involves intermediate particles of different flavors that can mix 
and oscillate. We discuss the implication of such oscillation, and in particular its interplay with the 
non-zero widths of the particles. We derive explicit formulae for differential decay rates involving 
both non-zero widths and oscillation, and show that in the case where the mass difference between 
the intermediate particle is of the order of their widths, both oscillation and width effects are 
important. An examination of the physical observables contained in these differential decay rates 
is provided. We calculate differential decay rates for cases in which the intermediate particles are 
either scalar s or fermions. 
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I. INTRODUCTION 



Many new physics scenarios predict cascade decays of new heavy degrees of freedom into 
Standard Model (SM) particles. A canonical example is the cascade decay of a squark into 
a quark plus two leptons and a neutralino q — > q\\ ill QllXv ft is well-known that for 
a cascade decay of the general form 

A -)■ XB -)■ XYC , (1) 

where X and Y are massless SM particles and ttia > itlb > mc, the differential decay rate, 
dTA/ds, possesses a kinematic edge located at |l|-|6( 



{m 2 A - m 2 B )[m 2 B — in 



c) 



ni 



S=(p X +py) 2 . (2) 



B 



This kinematic edge is in essence a step function in the differential decay rate distribu- 
tion, and it arises due to kinematic upper bounds on the on-shellness of the intermediate 
exchanged particle, B. The location of the kinematic edge provides an indirect means to 
either measure or constrain the masses of the A, B and C particles involved in the cascade. 
This mass measurement technique is called kinematic edge or endpoint method 0-13]. It 



is particularly important in the case that C is invisible, in which case the particle masses 
cannot be measured directly. 

In order to derive the kinematic edge in Eq. (j2J), one must assume that B is an on-shell 
mass eigenstate. This is a natural and plausible assumption to make, but it neglects the 
fact that B must also have a non-zero width, Tb > 0. One expects a non-zero width for 
B to smear the kinematic edge, because such a width smears out the invariant mass range 
within which B can be on-shell. However, for all phenomenologically important scenarios 
V b <C tub, so this smearing effect is considered to be small, and for this reason the role of Tb 
has been usually neglected. In some previous analyses, Tb has been incorporated into the 
differential decay rate by convolving the kinematic edge with a Breit-Wigner distribution 
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In many well-motivated theories the field B has not one but several flavors, which means 
that B is a superposition of multiple mass eigenstates that may mix together and oscillate. 
For example, this scenario is predicted in various SUSY theories 15|- and many proposals 



of ways to measure mass splittings, mixing and oscillation have been presented previously 
For just two flavors, it is well-known that the importance of the interference - 



the oscillation — between the mass eigenstates, denoted B\ and B 2 , is characterized by the 
dimensionless parameter 

x = ~y~ > ( 3 ) 

that is the ratio of the Bi 2 mass splitting to their average decay rate. In the case that i C 1 
or i > 1, oscillation is respectively unimportant because the oscillation length scale is too 
long or the oscillation is washed out. Oscillation effects, however, are significant in the case 
that x ~ 1. Due to the dependence of x on both the mass splitting and the decay rates, 
non-zero width and flavor oscillation effects cannot be independently considered. In other 
words, analysis of flavor oscillation requires the incorporation of the non-zero B widths into 
the computation of the differential decay rates. 

If the flavors do not oscillate significantly, or if interference is negligible due to x ^> 1, then 
we simply expect dT^/ds to feature multiple, distinct kinematic edges, each corresponding 



2 



to a single mass eigenstate, and the role of the widths should be unimportant. A detailed 
analysis of the physical information contained in the differential decay rates for the limit 



x — > oo, in which oscillation and widths are negligible, has been conducted in Ref. [25 
However, if oscillation is significant, then we not only expect interference terms to become 
important, but we also expect the form of dT^/ds to be smeared by the non-negligible 
widths. 

The purpose of this paper is to examine this prediction and its consequences in detail. We 
do this for the case of two-flavor mixing, with £? 12 either scalars or spin-1/2 fermions that 
interact with the external particles A, C, X, and Y via Yukawa-type interactions. In Sec. 
|II]we first redevelop the kinematic edge formalism via a field theoretic approach, accounting 
for the finite width with a Breit-Wigner propagator for B. For the simple case of a scalar c/> 3 
interaction with a single intermediate B, we show explicitly how the kinematic edge arises in 
the r# — > limit. We also show that the kinematic edges have their own well-defined 'edge 
width' that is a function of T B /mB, permitting us to quantify how much the kinematic edge 
is smeared for a given B width. 

In Sec. II III we introduce two-flavor mixing for both the scalar and fermion cases, and 
present the corresponding explicit results for cT^/ds in detail. As expected, we find that 
interference between the mass eigenstates is important only for the regime x ~ 1. We also 
rederive the result that in the fermion case, spin correlations between X and Y alter the 
shape of dTA/ds dramatically compared to the scalar B case, and verify that in the case 
of a vectorial coupling, the two-flavor fermionic case corresponds to the scalar one (see e.g. 
(26|). We present the explicit and detailed derivations of the results presented in Sec. IHIlin 
Appendix [A] 

Finally, in Sec. [IV] we briefly explore the physical observables contained in the differen- 
tial decay rates with oscillation, and how these may be used to constrain or measure the 
oscillation parameters, the masses, and the decay rates. We also propose a kinematic edge 
resolution criterion, which specifies under what conditions the two kinematic edges can be 
distinguished. We show that in most regions of parameter space the edge resolution crite- 
rion is simply x > 1, which aligns with our expectation that the edges should be resolvable 
when interference is negligible. We consider as an example in Sec. IIVI the special case that 
B consists of two flavors of sleptons, with oscillation parameters as motivated by gauge 
mediation SUSY breaking theories. Along with the kinematic edge constraints, we show 
that: the degree of oscillation — the magnitude of x — can be determined directly from the 
widths of the edges; the ratio of the kinematic edge step heights provides three observables, 
which uniquely constrain a parameter subspace involving x, the mass splitting and the mix- 
ing angle; and that the s = intercepts provide two other physical observables, which are 
measurable even if the edges cannot be resolved. Lastly, we show that for the alternative 
case that the Bs are fermions, the parameter space is constrained by these observables in an 
identical fashion to the scalars, the only difference being that the fermion parameter space 
is enlarged by one extra dimension compared to the scalar case. 



II. NON-ZERO WIDTH 

In this section, we examine the role of a non-zero width r#, and how it affects the sharp 
kinematic edge of Eq. ([2]). The usual derivation of Eq. (T5]) requires three assumptions: X 
and Y are massless; energy-momentum conservation; and, crucially, that B is on-shell. One 
may then derive Eq. (T5]) from kinematics alone. To include the finite width, we must instead 



3 




FIG. 1. Decay amplitudes for the particle A in the case that the intermediate states B\ 2 are either 
(a) scalars or (b) fermions. 

perform an explicit field-theoretic computation of the differential decay rate dTA/ds. 
A. Breit-Wigner Approximation 

In this paper, we will be mainly concerned with the differential decay rates associated 
with the two amplitudes shown in Fig. (TJ in which the intermediate exchanged particles Bi 2 
are respectively either scalars or fermions, and have Yukawa-type couplings to the external 
fields A, X, Y and C. We assume that X and Y are both leptons, which are denoted by 
the fields i a . We also assume, unless stated otherwise, that X and Y are massless 

m x = m Y = , (4) 

while the particles A, B 12 and C have cascade-ordered masses 

m A > mi )2 > m c ■ (5) 

To simplify notation we drop henceforth the subscript A from the decay rate, so that T A = T. 

Let us consider the propagator for B in the field-theoretic approach. Throughout the 
paper we shall always assume that 

T B /m B < 1 . (6) 

It is well-known that in the small width regime fl6]), and provided p B ~ m B , the exact 
propagator for an unstable particle is well-approximated by the Breit-Wigner propagator. 
(In a formal way, this approximation amounts to the leading order term of the propagator's 
Laurent series expanded around the pole p 2 = m 2 — imT.) In the same spirit as the sad- 
dle point approximation, integrals of amplitudes involving the propagator — precisely the 
objects with which we are concerned in this paper — are well-approximated by the Breit- 
Wigner function, provided the domain of phase space integration includes the Breit-Wigner 
maximum p 2 B = m 2 B . The reason is that the dominant contribution to such integrals is then 
due to an interval containing p\ = m 2 B , and the Breit-Wigner is the leading order approxi- 
mation within this interval. Outside this interval, the contribution to the integral by both 
the propagator and the Breit-Wigner are negligible, so we can approximate the propagator 
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by the Breit-Wigner for all p 2 B , even though the Breit-Wigner tails may be poor approxima- 
tions to the exact propagator tails. This argument extends naturally to objects involving 
sums of different propagators, that is, the differential decay rate for multiple flavors. 

Inclusion of the Breit-Wigner maximum p\ = m B within the domain of integration is 
guaranteed by assumption ([5]). Hence for any number of flavors, finite width effects in the 
differential decay rate are approximated to leading order in Tb/^ib by use of the Breit- 
Wigner propagators. 

For just two flavors, note that this result holds for all values of x. In particular, it 
is well-known that if the maxima are well-separated such that only the propagator tails 
overlap — i.e. x is large — then the Breit-Wigner approximation is a poor approximation 
to the corresponding interference term. However, we emphasize that if only the tails of 
the propagators overlap, then the corresponding interference term is also always negligible 
compared to the other terms contributing to the differential decay rate. As a result the failure 
of the Breit-Wigner approximation in the interference term for large x results in a negligible 
correction to the overall differential decay rate. So the Breit-Wigner approximation is a 
valid leading order approximation for all x. 

For the scalar, the Breit-Wigner approximation is 



Dsc(p 2 ) = Y-—. = , (7) 



while for a fermion 



f + y/m 2 — imY i> + m — iY/2 

D W = ' 1 2 2J _- v - ll 2 — rrr ' 8 

f p z — m z + imY p z — m z + imY 

due to Eq. (j^D- Note that we keep the iT/2 term as it is not necessarily small compared to 
ft + m, and in fact it produces leading order contributions in the case of intermediate particle 
oscillation. 



B. 3 Interaction 

To study the smearing of the kinematic edge due to T B ^ 0, it is instructive to first 
consider a toy 3 interaction calculation in which all the particles involved in the cascade 
are scalars and computational technicalities are therefore simplified. For further simplicity, 
we assume B is just a single scalar mass eigenstate. 

Consider such a 4> 3 scalar interaction, of the form 

£ = 9x4>a4>b4>x + gY4>c4>B4>Y ■ (9) 

The amplitude for the decay of A is 

lM = n> tyi z I im F • (10) 

Pb ~ m B + lm Bi- B 



Squaring this amplitude and integrating over phase space, one finds 

77=77+ (s) 

(11) 

7?=7?-0) 



dT_ gWv 1 tan -i^ 



ds 32(27r) 3 m 3 4 m B r B \T B 
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where 



r)±(s) = 1 + £(s) ± y£(s) 2 — m\m^/m 
5 



4 



We report the details of this calculation in Appendix |A] 

It is important to note that Tb oc <?y/ m s {gx,Y have mass dimension one here), so in 
the limit Tb/ttlb — > 0, g\jmB^B is finite. This means that Eq. (ITTj) has a well-defined 
normalization in the I" 'b/^b — > limit. To encode this explicitly, we therefore write 



2 



^ , (13) 



which is dimensionless. Adopting this convention, Eq. ( II ip becomes 

(14) 



^£ __9xd_, -i fm B \ 
ds 32(2vr)3m3 tan yTg J 



V=V-(s) 



which is manifestly finite in the limit Tb/tub — > 0. We will employ similar redefinitions of 
the coupling at the Y vertex throughout the paper. 

Overall momentum conservation in the A rest frame requires that s = m\ + m c — 2Ecrn? A) 
so it must always be that 

s < (m A - m c ) 2 = s max . (15) 

Observe that 77-(s max ) = ?7+(s max ), and r]±(s) become complex for s > s max , so the differential 
decay rate is always precisely zero for s > s max as expected. We emphasize that this 
maximum is not related to the kinematic edge, but is a distinct kinematic constraint. 

Let us now extract the kinematic edge from the differential decay rate. Expanding Eq. 
(JHJ) in rg/mg C 1 we have 



dT 7r 
~ds^2 



(sgn[r/ + (s)] -sgn[r?_(s)]} - ( — |— TTj— + ••• • (16) 



The leading order term produces a step function in s: For sgn[?7 + (s)] = — sgn[?y_(s)] the 
leading order term is 7r, while for sgn[^ + (s)] = sgn[?7_(s)], the leading order term is zero. 
The location of this step, Sq, must therefore satisfy either 

77 + (s ) = 0, or 77_(s )=0. (17) 

(The ratio of to the geometric mean of and mc determines which function is zero 
in Eqs. f|T7|) . If m^mc /m 2 B < 1 then ?7_(s ) = and r] + (s ) > 0; if mA^nc / m B > 1 then 
V+( s o) — and ^-(s ) < 0; while if m A mc /m 2 B = 1 then both rj±(s ) = 0.) 

Using the definitions ffl2l . one may verify that the solution to Eqs. (fTT|) is always 

_ (m\-m B )(m B -m 2 c ) 

S — 2 ' \ l °> 

m B 

which is precisely the expected kinematic edge fl2]). We have therefore shown that the zeroth 
order contribution in Yb/ttlb to Eq. f fl6|) produces the kinematic edge, while terms of the 
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FIG. 2. Differential decay rate with parameter choice m^/m^ = 2, mc/rriB = 0.2 and Tb/ttlb = 
1CT 1 (solid) or Vs/ms = 10~ 3 (dashed). The kinematic edge so/m B = 2.88 clearly emerges as 
T B /m B ->■ 0. 



order T B /m B and higher smear the edge into Eq. (JHJ). A plot of the differential decay 
rate for different r^/m^ is shown in Fig. [20 Notice that if m 2 B = rriATnc, then s = s max , 
so the kinematic edge collides with the overall kinematic constraint. For all other cases 
so < s max . Note also that the terms of the order r b/^b and higher in Eq. (fl6l) each diverge 
at the kinematic edge, while the resummed expression in Eq. (TL41) is finite. This occurs 
because the formal expansion of the Breit-Wigner propagator in powers of T B / itib does not 
converge. Hence, even though expanding Eq. (|14p provides us insight into the kinematic edge 
structure, nonetheless the full expression ( I14jl is itself the leading order differential decay rate 
in Tb/ttlb- Similarly, throughout this paper we will not expand the closed form functions 
produced by the Breit-Wigner approximation, although we will expand their prefactors. 

That the kinematic edge arises from the T B /m B — > limit of an arctangent function is 
the main result of this section. We will see below that these functions, and their associated 
kinematic edges, are a general feature of the differential decay rates considered in this paper. 



C. On-shellness of B 



It is also instructive at this point to consider how the on-shellness of B is encoded in this 
field-theoretic derivation of the kinematic edge. The amplitude squared can be written as 

1 1 



\M\' 



9x9y 



9x9y 



(p 2 -pI)(p 2 



vf 



_ r _ p 2 p 2 2 ' 

indicates complex conjugation. The key observation is 
that if one takes the Tb/ttlb — > limit before rather than after computation of dT/ds, one 
finds that the amplitude squared is 



where j>\ = m B — ims^B and the 



hm 



Im 



f 



p z — m 



B 



B 



2tt 



m 



9x9yKv ~ m B ) . 



(20) 



B 



1 Readers expert in the kinematic edge method may wonder why the differential decay rate in Fig. [2] is 
rectangular in shape, rather than the usual triangle. The reason is that we have plotted here dT/ds rather 
than dT/d^/s, the latter convention being common in the Literature because the background is usually 
flat in -y/s. However, the former convention is also used (see e.g. 
always consider dT/ds. 
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26l|)- Throughout this paper we shall 



We then see that the differential decay rate in the limit Tb/tub — s-Ois simply a phase space 
integral of a delta function that forces B to be on-shell. This integral yields a step function 
as expected, with the edge occurring at the value of s for which B can no longer be on-shell. 



D. Edge Width 

Let us now specify how the smearing of the kinematic edge is related to non-zero Tb/tub- 
In order to characterize the amount of smearing of the kinematic edge, one may show that 
the gradient of the differential decay rate, d 2 T/ds 2 , is of the form 

d^/ds 2 = /( g - , (21) 
(s - s ) 2 + a 2 

where f(s) is a smooth function of s that is slowly varying compared to the Breit-Wigner 
factor near the kinematic edge. Hence near the edge 



d 2 T/ds 2 * , ff*jL - , (22) 



in the same spirit as the saddle point approximation. This Breit-Wigner is clearly maximal 
at the kinematic edge and has full width at half-maximum (FWHM) 



a ~ 2m B 



2 



m B 



(T B /m B ) 2 + 1 - m 2 A m 2 c /m B . (23) 



We call a the edge width. Note that in general (TB/mB) 2 is not necessarily small compared 
to 1 — m 2 A m 2 c jm B - For the examples shown in Fig. [2], we have cr/m 2 B ~ 0.2 or 2 x 10 -3 
respectively, which match the naively expected orders of magnitude. The FWHM of d 2 T/ds 2 
is in principle a measurable quantity: Measurement of a provides constraints on the size of 
the ratio Vs/mB, rendering Eq. ( 1231) as an important results. 



E. Non-zero Width of A 

So far in this discussion we have treated A as an on-shell external state. In practice, 
however, A itself is an intermediate state of a yet larger cascade that started with a heavier 
mother particle, denoted A'. In such a scenario, one measures dT A i/ds instead of dT A /ds, 
and since A has a non-zero width, then the amplitude corresponding to this differential 
decay rate has the form 

( ! ^ ( ! ^ . (24) 

\p 2 A -m 2 A + %m A Y A ) \ (va ~ Vx) 2 - m% + im B T B ) 

Here the smearing due to T A convolves with that of V b'- If is sufficiently large then we 
expect the B kinematic edge structure to be lost — i.e. smeared away — by this convolution. 
For non-zero T A , not only does this convolution generally further smear the kinematic edge, 
but it also means we cannot generally distinguish the smearing effects of T A from in the 
differential decay rate. That is, the results of Sec. Ill B\ which account for smearing due to 
non-zero T B alone, would be invalid. 
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Our analysis of non-zero T b effects in this section implies that smearing due to Ta vanishes 
as the parameter T^/m^ — > 0. In particular, from the amplitude ()2l]) we expect the Ta 
smearing to be negligible if 

T A /m A < T B /m B . (25) 

For multiple flavors we require V a/^a *C ^j/ m j f° r an 3- 111 this regime, the effect of A's 
non-zero width on each kinematic edge is negligible compared to the effects of the respective 
non-zero B widths. 

From Eq. (THj) . we have 

T A g\~g\ s < g\g\ , , 

mA rn A rn A m A 

since s /m A < 1. So smearing due to A's non-zero width is negligible provided the coupling 
gx is sufficiently small. For the remainder of this paper, we shall always assume Eq. (T2"5"j) is 
satisfied for all flavors, so that smearing due to A is negligible. 

We shall now use the insight we have gained into finite width effects from this simple 3 
theory to study the kinematic edge with flavor oscillation. 



III. FLAVOR OSCILLATION 

In this section we present results for the cases that B\ 2 are scalars or fermions, which we 
call the intermediate scalar and intermediate fermion cases respectively. 



A. Oscillation Parameters and Small Width Regime 

Before proceeding, let us define the following usual oscillation parameters in terms of the 
Bi t 2 mass and decay rates, and Ti^: 

_ m 2 + mi _ - _ Ta + Ti Ar _ r r {0 v\ 

m = , Am = m 2 -mi , T = , Ar = T 2 - T 1 , (27) 

where m 2 > m% and 

Am Ar Am . . 

x = , y =^ z = . (28) 

r ' y 2T ' 2m v ; 

Note that x > is unbounded, while — 1 < y < 1 and < z < 1. The small width regimes 
for each mass eigenstate are defined, as usual, by Tj/rrij < 1. It is convenient to define the 
parameters 

_ -r / _ r / mid + m 2 e 2 

e,; = I ,7m, , e = I /m = . 

3 31 3 ' ' m 1 +m 2 

Just as before, we shall always assume small widths tj 1. Observe that since m 12 and 
ei j2 are positive definite quantities, then this assumption implies e < 1. 

The four parameters m, x, y and z are independent so they uniquely specify m^, Ti i2 
and ei j2 , viz. 

2mz , . 2z 2zfl±y\ 

m 1>2 = m(l±z), r lj2 = (l±y , e = — , e lj2 = — — . 29 

X X X \ l±z J 



9 



It is clear that in the small width approximation 

z < x. (30) 

Ideally, we may present all the differential decay rates just in terms of the oscillation pa- 
rameters m, x, y and z alone. However, for the sake of compactness and clarity, we shall 
instead present our results in terms of a mixture of both m, x, y and z as well as m^, Ti^, 
e and £1,2 with the understanding that the latter may be expressed in terms of the former 
via Eqs.' flZ9) . 

Finally, we assume CP conservation. Thus, for two flavors the mixing matrix U, as 
defined in what follows (cf. Eqs. ( 152|) and (|42p ). is real orthogonal and has a single physical 
mixing angle, 9. We write 

U= 1 ( C ° S ° S[n °) ■ (31) 
2 V — sin 6 cos V J 



B. Intermediate Scalars 

First consider the differential decay rate due to the amplitude in Fig. [Th. Here B is a 
superposition of two mass eigenstates with B\ and B2 scalars of two different flavors, while 
A, X, Y and C are fermions, with Yukawa-type vertices defined by 

C s = $ A (g?P L + glP R )t a U^^ B + ^c{glP L + g Y R P R y a U ai ^ B ■ (32) 

We will not report here the study of a single intermediate scalar for this interaction. As 
shown in Appendix this case does not differ considerably from the 3 case studied in 
details in Section III Bl The only difference is the presence of linear and logarithmic terms 
in s, which are suppressed by factors of e. 

We write the differential decay rate as 




The first two terms come respectively from the squared single B\ and B 2 contributions. The 
final term is the interference term for the two different flavors. The subscript 'sc' denotes 
an internal scalar. Note also that just as for the 3 case in Eq. ([13]) . the combination 
(g\ R ) 2 m/r is finite in the zero width limit (the couplings are now dimensionless) , so we 
write, 

{Slnfj = (9ln) 2 , (34) 

in which g^ R are finite. In Eq. fl33|) we have already removed the 1/e = m/T factor absorbed 
by g\ R and replaced it with g\ R . 

To leading order in e the square terms are 

ri=rf + 

(35) 



ds 



\U a i\ 2 \U^\ 2 s J fl + 



(27T) 



m" 



tan 
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in which we have defined (cf. Eqs. ( IT2|) and ({18]) ) 



m 2 )(m 2 



m, 



c) 



Us) = 



m 2 

s — (m 2 A + m 



c) 



2m] 



m\m 2 c /m 4 : 



(36) 



but 



In Eq. (|35l) we have discarded terms whose coefficients are subleading order in e or e^, 
just as for the 3 example we have not expanded the arctangent function itself, in order to 
avoid creating artificial divergences at the kinematic edges. 
The interference term is 



ds 



m 



(27T) 



U al U^U a2 U^ 2 
x 2 + 1 



^ j^etan- 1 


v 


+ £f c log 


2 i 2 


- 1 j=l,2 I 









ri=rn 



V=V1 



(37) 



with coefficients that are to leading order given by 

(m\ — m 2 )(m 2 — tuq) 



BL 



rrr 

-IV -AO 

l > 2 sc 



(38) 



Several comments are necessary concerning the interference term and coefficients pre- 
sented in Eqs. ( 1371) and ( f38l) . We emphasize first that the coefficients are valid only to 
leading order in e. We have also dropped contributions of 0{z) and higher to the interfer- 
ence term prefactors, because their contributions are always suppressed. To see this, note 
that the interference terms can be written generally in the form 



dY 



a/3 
12 



ds 



a + bz 



(39) 



where a and b are arbitrary linear combinations of the arctangent and logarithm functions 
and a/b ~ 1. It is clear that the denominator ensures the interference term is relevant 
only for x < 1: If instead then the denominator suppresses the entire interference 

term, so the O(z) terms are certainly unimportant. In the case that x < 1, then z C 1 
by Eq. (|30|) . Hence the O(z) terms are always negligible. We emphasize, however, that 
we have not discarded the z dependence of the arctangent and logarithm arguments, only 
that of their prefactors. The ability to drop the O(z) contributions in the interference term 
prefactors applies similarly to the intermediate fermion case. The full expressions, including 
O(z) terms are reported in Appendix IA 21 

The square terms in Eq. fl35|) bear obvious similarities to the 3 result in Eq. (j!4p . The 
same analysis applies. Each square term becomes a step function in the zero width limit, 
producing a kinematic edge at 



im. 



m 2 )(m 2 



m 



c) 



(40) 
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FIG. 3. Left: Intermediate scalar differential decay rate for a = j3 = 1, with the parameter choices 
6 = 7r/4, niA/m = 2, mc/m = 0.2, z = y = 0.1, and x = 4 (solid) or x = 4 x 10 2 (dashed). 
These correspond to e = 5 x 10~ 2 , 5 x 10 -4 respectively. The kinematic edges at Sq/to 2 = 2.7 and 
s^/m 2 = 3.0 are smeared as e becomes larger. Right: Intermediate scalar differential decay rate for 
the same parameter choices and x = 4 with (solid) and without (dashed) the interference terms. 



respectively, but these edges are smeared out by non-zero 6j. Eq. (|2"§|) implies that for a 
fixed z and y, the smaller x is, the larger the Sj become. Hence in the oscillation regime, 
x ~ 1, the kinematic edges are more smeared compared to the x ^> 1 case. Similarly, the 
arctangent contributions to the interference term have edges which are smeared according 
to the size of the 6j parameters, with edge widths given by Eq. (1231) . As expected, strong 
flavor oscillation smears the kinematic edges. An example of this smearing of the kinematic 
edges is shown in Fig. [3] for fixed z but varying x. In this figure we have chosen z = 0.1, a 
relatively large value, in order that the edges are visually distinct. 

As mentioned above, the interference term has a 1/ (1+x 2 ) prefactor, which means that for 
x 3> 1 the entire interference term is suppressed. Since x ^> 1 corresponds to no interference 
due to flavor oscillation 'wash-out', this is precisely the expected behavior. In Fig. [3] an 
example of the relative importance of the interference term contribution is shown graphically 
for the case that x ~ 1. As can be seen, for the chosen parameters the interference terms 
appear to slightly enhance the sharpness of the edges. Obviously, as x becomes smaller, the 
role of the interference term becomes more significant. 

Another new feature of the interference term are the logarithm terms. These logs do not 
contain kinematic edges. Instead, in terms of e for fixed z, and being careful to include the 
l/(x 2 + 1) and B sc ~ x prefactors, the log terms are of form 



dT 



a/5 
12 



ds 



sc,log 



ze 
e 2 + z 2 



log {rf + e 2 ) 



(41) 



From Eq. ( fTTl) . at a kinematic edge we have either n± = 0, so the log terms are largest at a 
kinematic edge. However, in the e — > (or x — > oo) limit these log terms manifestly vanish. 
Notice also that for z — > with x fixed, the large logs at the kinematic edge cancel due to 
the different signs of B J SC . 
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C. Single Intermediate Fermion 



In contrast to the intermediate scalar case, a single intermediate fermion case presents 
novelties which are worth discussing. In particular, since the B is now a spin- 1/2 particle, 
we expect spin correlations between X and Y . These correlations lead to a differential decay 
rate with triangular or trapezoidal envelope rather than the square envelope found for the 
intermediate scalar case. 

Consider the following interactions for a single intermediate fermionic B, 



(42) 



Squaring the amplitude in Fig. [Tb, summing (averaging) over final (initial) spins, and 
integrating over phase space, one finds differential decay rate of the form 



dT 

ds 



[< 



«M) 2 



dT 



+ [(^) 2 + (^D 2 ]^ 



(43) 



The subscript 'f denotes an internal fermion. The couplings g\ R are defined as in Eqs. ( TT3]) 
and (134"]) . and the ttib/Tb factor has already been similarly factored out of the decay rates 
dT/ds\±. The two terms in Eq. (143]) arise from two purely chiral interactions: i.e. the cases 
9l = 9r = ^ oi 9r = 9r = respectively. We therefore call dT/ds\± the chiral differential 
decay rates. The chiral differential decay rates turn out to have respectively a positive or 
negative slope in s, whence the subscript (see Appendix [A] for details). 
One finds to leading order in T B /m B <C 1, 



dT 

ds 

dT 

ds 



s 



(2vr) 3 m3 



(27r) 3 m^ 



tan 1 



tan 



m B 



V 



m B 



r)=r) + (s) 

T)=T)-(s) 

V=V+(s) 
ri=ri-(s) 



(44) 



where r]± and so are defined in Eqs. ( fT2l) and f|T8|) respectively. 

A few remarks about Eqs. (j43|) and (1441) are in order. First, obviously the linear s 
dependence of the arctangent prefactors in Eqs. (PHI) changes the overall shape of the chiral 
differential decay rates from a rectangle into a triangle of either positive or negative slope. 
Second, just as for the intermediate scalar case, the arctangent functions produce a kinematic 
edge at s = So which is smeared by non-zero r#. Note, however, that the negatively sloped 
differential decay rate is precisely zero at the kinematic edge. Finally, from Eqs. (14"4"|) and 
( l4"3"j) . we see that in dT/ds\{ the linear s dependent piece of the arctangent prefactor is 
proportional to 

M~9l? + i9 X R ~9l) 2 ] - Mfn? + {9 X R m ■ (45) 

This s dependence disappears if either the X or Y vertex has a vectorial coupling, i.e g£ = 
or g\ = gp[ respectively. In this case we obtain a square envelope for the differential decay 
rate. In other words, if at least one coupling is vectorial then there are no spin correlations 
between X and Y, provided one averages over their initial and final spins. This means 
that for such vectorial coupling, one cannot distinguish a single intermediate fermion from 
an intermediate scalar using the shape of the differential decay rate. This is a well-known 
result (see e.g. 26]). Plots of the chiral differential decay rates are shown in Fig. HI 
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FIG. 4. Chiral differential decay rates for the cases = g\ = (left) and g^ = g\ = (right) 
with parameter choice mi/m^ = 2, mc/ms = 0.2 and T^/ms = 1CP 1 (solid) or Vs/niB = 1CT 3 
(dashed). For the positive slope distribution, the kinematic edge so/m 2 B = 2.88 clearly emerges as 
r^/rriB — > 0, while the negative slope is always precisely zero at the kinematic edge. 



D. Intermediate Fermion with Flavor Oscillation 



Let us next consider the differential decay rate for the amplitude in Fig. [TJd with two- 
flavor mixing. This time there are two intermediate fermions Bi and B 2 , with Yukawa-type 
vertices 

£ f = $ B (gfP L + g£P R )U m *£ a ^ A + ^ B {g Y L P L + g Y R P R )U m *£ a ^ c . (46) 
Just as for the single intermediate fermion, the differential decay rate is of the form 



ds 



Mf L f + {9 X R ~9l) 2 ] 



dT a ? 



ds 



Mfnf + (<?K) 2 ] 



dT ali 



ds 



(47) 



Again the subscript 'f denotes an internal fermion. The chiral differential decay rates, 
dT/ds\±, now have both square and interference terms, so we write 



dT c 



ds 



dTl 



a8 



± 



ds 



+ 



dTf 



ds 



+ 



± 



ds 



± 



The square terms are similar to Eqs. ( Hi]) . To leading order in e, they are 



(48) 



dTf 



ds 



dTf 



ds 



ljjaj\2ijj/3ji2 
(2^)3 

\U aj \ 2 \U l3j \ 2 

(2^p 



m , 



l+(-)>z\ _! 

tan 



4 - s 



m 



A J 



i + {-yy 
i + (-y y 



tan 



v=v+ 



r]=rf. 



(49) 



The interference terms are similar to Eq. ()3 



a/3 
12 



(27r) 3 m^ 



U al U^U a2 U^ 2 
x 2 + 1 



<^ ^ tan" 1 

J'=l,2 



+ #jjog 



2 , 2 

77 



V=V- 



(50) 
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except that the coefficients (to leading order in e) are 





s 

9 ' 






= -n j 


X s 
2 m? ' 


aL 


— A? - 


s 

m z 


Bl 


= Bl + 


( ) jX S 
1 J 2m 2 



(51) 

We have again dropped the 0{z) terms, as discussed in Sec. IIIIBl Notice that the square 
terms dTj/ds\- are precisely zero at their kinematic edges and thereafter. As a result, rather 
than a double edge, the dT /ds\ - rate has two kinks, which are smeared as e increases. Plots 
of the chiral differential decay rates dT /ds\± and summed differential decay rates (f43|) are 
presented in Fig. [5j 

Our analysis in Sec. IIII Bl concerning the role of the logarithms, arctangents and l/(x 2 + l) 
factors in the interference terms for scalars applies equally to our results for the intermediate 
fermions. As mentioned above, the main difference between the intermediate scalar and 
intermediate fermion results are the linear s prefactors, along with the fact that there are 
twice as many terms corresponding to the two chiral coupling cases. Finally, note that 
in the case a coupling is vectorial, both the square and interference terms reduce to the 
intermediate scalar case: The linear s dependence in the prefactors manifestly cancels and 
the fermion and scalar couplings coincide, just as for a single mass eigenstate. 

This cancellation does not, however, persist in the log interference terms at higher e order 
(see Appendix I A 31 for details). This non-cancellation may be an artifact of the Breit-Wigner 
approximation [27, [28[ and is certainly negligible here. Yet, it is worth pointing out that 
such a difference could open possibilities for intermediate particle spin determination in the 
vectorial coupling case. A possible physical origin of the difference between the bosonic 
and the fermionic cases may be off-shellness effects in the interference term: When the 
intermediate particles are far off-shell, e.g. in the 7711,2 ^> mA > mc regime, the differential 
decay rate is very different depending on the spin of the intermediate particle. Since off- 
shellness is parametrized by e, to compute such off-shell spin effects would require analysis 
of higher order e terms. This is beyond the scope of this paper, and thus we do not exploit 
this possible avenue for spin determination any further here. 



IV. OBSERVABLES 

Having presented the differential decay rates for the intermediate fermion and scalar cases 
with flavor oscillation in Sec. IHH we now examine what physical observables can in principle 
be measured from these distributions. We emphasize that we do not study the feasibility 
of such measurements. A detailed analysis of the physical information contained in the 
differential decay rates for the limit x — > oo, in which oscillation and widths are negligible, 



has been conducted in Ref. 25 



For the purposes of this discussion, in this section we will focus on the canonical example 
of cascade decay with an intermediate slepton (that is, a scalar) 

X? « « *S • (52) 



15 



1.0 2.0 2.73.0 4.0 

s/m 2 



1.0 2.0 2.73.0 4.0 

s/m 2 




FIG. 5. Internal fermion chiral differential decay rates dT u /ds\+ (upper left), dT 11 /ds\- (upper 
right), and summed differential decay rates dT u /ds (lower left and right). The parameter choices 
are 9 = 7r/4, rriA/m = 2, mc/m = 0.2, z = y = 0.1, and x = 4 (solid) or x = 4 x 10 2 (dashed), 
which correspond to e = 5 x 10~ 2 and 5 x 10 respectively. For the summed rates we choose 
couplings such that (g* g\f + (g% g\) 2 = 1/2, (g* g^) 2 + (g%gl) 2 = 3/2 (lower left) and vice- 
versa (lower right). The kinematic edges in dV/ds\ + at s\/m 2 = 2.7 and Sq/to 2 = 3.0 are evident as 
x becomes large, while in dT/ds\- we see instead two kinks. The trapezoidal shape and kinematic 
edges of the summed rates dT 11 /ds are manifest. 



In line with our discussion so far, we shall assume that there are only two flavors of slepton 
which mix significantly. We shall also assume they couple to the neutralino Xi (X2) an d the 
electron e (muon /1) via CP preserving interactions of the form in Eq. (132|) . We adopt the 
notation 

?Pa,c = Xi,2 , i 1 ' 2 = e,/i . (53) 

In contrast to the notation denned in Eq. (|31|) . for the sake of clarity in this section we will 
henceforth assign the flavor indices of the mixing matrix to be a = e, [i. 

The SUSY spectrum depends on the mechanism of SUSY breaking. If the breaking is 
mediated in a flavor-universal manner then the mass splitting is small. For example, gauge 
mediation in SUSY breaking theories naturally gives mass splittings and decay widths for 
sleptons of the order GeV. Since existing bounds on the slepton masses imply m- t > 10 2 
GeV, then for this decay we consider oscillation parameters 

m ~ 10 2 GeV , y ~ , z ~ 10~ 2 , x ~ 1 . (54) 
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These parameters satisfy the small width condition <C 1, and flavor oscillation is strong. 

A. Parameter Counting: Kinematic Edges 

As explained in Sec. IIII A[ measurement of the four oscillation parameters m, x, y and z 
uniquely determines the slepton masses and widths. The other physical parameters we wish 
to measure are the mixing angle 9, the neutralino masses m xi and m X2 , and the couplings 
g^n and g\ R , so we have eleven physical parameters of interest in total. However, from Eq. 
f )33|) only the combination 

9 - Mf + {9 X R f] Wl? + (fn) 2 ] , (55) 

appears in the differential decay rate, so only g can be measured. Thus we have an eight 
dimensional parameter space. 

In the case of two-flavor mixing and with CP conservation, the differential decay rates 
are symmetric in flavor indices, as can be seen from Eqs. (1351) and ([57)1 ■ That is, dT al3 /ds = 
dT l3a /ds. As a result, there are three independent differential decay rate distributions which 
can be measured, namely 

dT ee dT en ^ 

ds ' ds ds ' ds 
Note that in the most general case U can be complex, and the interference prefactors then 
generally contain jj m u l3t *U aj *U l3j . As a result dT a/3 /ds ^ dT^/ds, and then each final state 
provides an independent differential decay rate. 

The positions of the kinematic edges ( H0)l are obviously lepton flavor independent, and 
therefore are the same for these three distributions. The edges therefore yield two indepen- 
dent constraints on the four- dimensional parameter subspace {m Xl ,m X2 ,m, z}, constraining 
it to a two-dimensional surface. The full parameter space is constrained by the edges to 
a six- dimensional surface. We now must seek other observables to further constrain the 
parameter space. 

B. Direct Measurement of the Widths and Oscillation 

In Sec. Ill Dl we established that the smearing of the edges is characterized by the edge 
width a, defined in Eq. ([23)) . The edge width approximates the full width at half maximum of 
the derivative of the differential decay rate, dT 2 /ds 2 , and is therefore a measurable quantity. 
Assuming that 1 — m 2 2 /m 4 ~ and since y ~ 0, z ~ 10 -2 <C 1, then for the 

intermediate slepton cascade the edge width reduces simply to 

a „ 2c m4 ~ « _ 4m 2 ^ / m 2 x m 2 X2 \ ^ 
m? x \ m 4 / ' 

for both kinematic edges. Since this is lepton flavor independent, all three differential decay 
rates feature two kinematic edges with the same edge width, yielding one independent further 
constraint on the parameter space. 

Significantly, since the magnitudes of m and z are fixed by other theoretical considerations 
(as in Eq. (jMJ)), and provided 1 —m xi m 2 /m l ~ 0(1), then measurement of a provides the 
magnitude of x. As a consequence, the measurement of the smearing of the edges — the 
edge width — measures the degree of flavor oscillation. 
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1.0 2.0 2.73.0 4.0 

s/m 2 

FIG. 6. Intermediate scalar differential decay rate, with the parameter choices 6 = tt/4, mA/m = 2, 
mc/m = 0.2, y = 0, z = 0.1, and x = 1 (solid), x = 2 (dashed) or x = 3 (dotted). These correspond 
to e = 0.2, e = 0.1 and e = 0.02 respectively. 



C. Edge Resolution Criterion 

So far in this discussion we have assumed that it is possible to distinguish the two edges, 
that is to say, it is possible to resolve them. But if the separation of the edges is of similar or 
smaller size than their widths, 01,2, then it is reasonable to assume that one may not resolve 
the two edges. This leads to a natural edge resolution criterion. For two kinematic edges to 
be resolvable, we require the separation of the edges to be greater than their average width, 
that is 

\4-4\>^- (58) 

For z < 1 and 1 — mL x mL „/m A ~ 0(1), this criterion reduces to just a simple restriction on 



XI X2 



x > 1 . (59) 



This is an interesting result, because the resolvability depends on both the separation and 
the widths of the edges, so naively we would expect the resolution criterion to depend on 
both x and z. That said, this result implies that oscillation and edge resolvability are 
inversely correlated, which aligns with our expectation that the edges should be resolvable 
when interference is negligible. For the intermediate slepton cascade we have x ~ 1, so 
the edges may not always be resolvable, but they are if we restrict our attention to x > 1, 
which is often the regime of interest. A graphical demonstration of the sensitivity of this 
criterion is provided in Fig. |6j For x = 2 the edges are already visibly resolvable. Note 
that the parameters have been chosen in this figure to match our prior discussion, while still 
satisfying 2,e C 1 and 1 — m^m^/m 4 ~ 0(1), and do not arise from the slepton cascade. 



D. Edge Degeneracy and the Geometric Mean 

Before continuing, we wish to point out that kinematic edges can be irresolvable even in 
the absence of strong oscillation. That is, there are special points in parameter space at 
which two kinematic edges fHOj) are degenerate, i.e. = Sq. This occurs not only for the 
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obvious case that mi = m 2 , but also if 

mim 2 = mArnc , (60) 

that is, if the geometric means of m 12 and m^c are the same. In terms of the parameter z, 
this condition corresponds to 

, = Ji-masc. (6i) 

V vnr 

If we require z < 1, then it follows that 1 — m 2 m 2 lm A <C 1, too. So Eq. (|60|) is satisfied 
in a different regime of parameter space to that considered in Sec. IIV CI 

The physical origin of this degeneracy can be understood by observing that if X and Y 
are massless, then the invariant mass s = 2(\px\\pr \ — Px-Py) < 4|px||pV|- One derives the 
kinematic edge, s J , from kinematics by maximizing s subject to Bj being on-shell. It is clear 
that a given value of s can be obtained from many configurations of px and py. For instance, 
s has the symmetry \p x \ — > Mpy\, \py\ \px\/K VA G R. It is therefore not surprising that 
there exist different on-shell B masses with correspondingly different px,Y configurations, but 
with so the same. In particular, it can be shown that mi — > mAvncl m \ — m 2 is equivalent 
to A = m^/mc. So if mim 2 = m^mc, then the two kinematic edges are degenerate. 

When Eq. (160]) is satisfied, the kinematic edges collide, so that one may not determine 
the number of intermediate mass eigenstates from the differential decay rate dT/ds. There 
may, of course, still be oscillation between the two B's in this case or none at all. The moral 
of this section and Sec. IIV CI is that the observation of a certain number of kinematic edges 
provides only a lower bound on the number of intermediate degrees of freedom. 

E. Step Height Ratios 

If the edges are resolvable, the ratio of the heights of the step functions that form the 
kinematic edges is another observable. For y ~ and z < 1, using the results of Sec. IIII Bl 
one may show that this step height ratio is to leading order in z 



R af3 = 



\u al 


2ppi\2 + U al U pi U a2 U P2[ 1 + 


(x/ir) log 


■e}/(x 2 


+ 1) 


pea 


2|[//32|2 _|_ JJalJjpiUcau P2^l _ 


(x/tt) log 


;e]/(x 2 


+ 1) 



(62) 



The step height ratio is lepton flavor dependent, yielding three constraints on the three 
dimensional {z, x, 9} parameter space. Consequently, the step height ratios probe z, x and 
9. Reversing the parameter counting above, measurement of a now constrains m uniquely, 
while the two kinematic edges are left to constrain m Xl and m X2 . With y ~ 0, only the 
coupling combination g is left unconstrained by these observables. 

F. s = Intercepts 

We now consider the value of the decay rate at s — 0, which we call the s = inter- 
cept. These intercepts are another three physical observables. Since they depend on overall 
normalization, they permit measurement of g, and moreover, they are measurable even if 
the edges cannot be resolved. Explicitly, for the three independent differential decay rates, 
these intercepts are to leading order 

T ap = ~ ( m xi ~ m 2 ){m 2 - m 2 X2 ) 
8n 2 m A m 2 



pal\2pPl\2 + pcx2\2pP2f 



X 2 + 1 



-U^U^U" 2 ^ 2 



■ (63) 
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Unlike the step height ratios, not only does I e>1 = J Me , but also I ee = J MM , so the inter- 
cepts actually provide just two independent observables. The kinematic edges, edge width, 
step height ratios and s = intercepts thus provide eight physical observables for the 
original eight-dimensional parameter space. Alternatively, if we set y ~ 0, then the seven- 
dimensional parameter space is over-constrained. Note also that the single independent 
ratio of these intercepts is normalization independent, and therefore constrains just the 
{x, 9} parameter subspace. 

G. Intermediate Fermion Observables 

Let us now consider the physical observables for the case of two-flavor mixing with inter- 
mediate fermions. It is convenient to define 

~9+ - {9 X l~9 Y l? + {9 X R ~9l? , ~9- - {9 X L ~9l? + (gMT , (64) 

which are the measurable coupling combinations analogous to g above. The parameter space 
for fermions is therefore enlarged by one dimension compared to the scalar case. We assume 
the same oscillation parameter structure as in Eqs. (|54|) . 

Just as for the intermediate scalar case above, there are still three independent differential 
decay distributions (j5B"|) ; the kinematic edges and edge width still provide respectively two 
and one physical observables that constrain the parameter space; and the edge resolution 
criterion (loTJj) is the same. Apart from the dimensionalities of their parameter spaces, the 
only other difference between the scalar and fermionic cases is that for the latter, rather 
than step height ratios one must instead measure the ratio of trapezoid apex heights. That 
is, one must extrapolate the smeared edges into sharp corners — the solid to the dotted 
lines in Fig. [5] — and measure the ratio of the trapezoid apex heights so formed. This is 
obviously possible, but will presumably introduce more error into the experimental results. 
At leading order in e and z these apex height ratios are identical to the step height ratios 
in Eq. flS2D, and the leading order contribution is due to dT/ds\ + only. 

The s = intercepts are similarly the same as in Eq. (153]) but with g replaced by just 
<7_, since dT/ds\ + is manifestly zero at s = 0, so the intercepts depend only on g_. One 
can constrain g + by measuring the apex heights themselves, rather than their ratio, which 
depend only on g + to leading order in z. 

V. CONCLUSION 

In this paper we have presented explicit differential decay rates for cascade decays of 
the form ([I]) as a function of the XY invariant mass, s, including the effects of both finite 
particle widths and flavor oscillation. In particular, we considered both scalar and spin- 1/2 
intermediate particles with Yukawa-type vertices. Our results successfully reproduce the 
usual kinematic edge results in zero width, zero oscillation limit, and we have shown how to 
quantify the degree of kinematic edge smearing for finite widths. 

The main results of this paper, however, involve the analysis of the interplay between 
non-zero width and flavor oscillation effects, as characterized by the parameters x = Am/f 
and z = Am /2m. Not only does x control the degree of oscillation or interference between 
the two B mass eigenstates, it also plays an important role in the degree of smearing of 



20 



the two kinematic edges and their resolvability. In general, the larger x, the smaller the 
smearing and interference, and the greater the resolvability. 

We have also shown how various physical observables of the differential decay rates can 
be used to constrain the oscillation parameter space, mainly for the special case of a slepton 
cascade with parameters motivated by gauge mediation SUSY breaking theories. In partic- 
ular, apart from the kinematic edges, these observables include the edge widths, step height 
ratios and the s = intercepts. Building on the special case considered in Sec. HV] a subject 
of future work may be to apply our explicit results to more realistic, more complicated or 
more general scenarios. 

Another avenue of study may be to extend the current two-flavor treatment to the three 
or more flavors. Furthermore, one might lift the CP conservation assumption which we 
have made throughout this paper. Such CP violation will produce CP violating interference 
terms, which may have interesting physical effects. We plan to address these issues in the 
future. 
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Appendix A: Phase Space Integral and Non-Zero Width 

In the following we compute the differential decay rate dT/ds for a cascade decay of the 
general form ([T]) 

A -> XB -> XYC . (Al) 

Although some of what is reported in this appendix is already well-known, it is our hope 
that its recapitulation here is a useful reference to the reader. 



For the three body decay (1A1J the differential decay rate is |29|, |30 



r-p \M\ 2 d 3 p x d 3 p Y d 3 p c (4) 

dT = ^]^~ A 2E^2E^2E^ 6 {PA -Pi-*- Pc) , (A2) 

where Ej = (p| + m 2 ) 1 / 2 and \Ai\ 2 is the amplitude-squared, to be specified later. The 
differential decay rate is obviously defined in the A rest frame. The amplitude-squared \ Ai\ 2 
is averaged over all spin indices (if any). As a result, note that \Ai\ 2 has spherical symmetry 
under arbitrary spatial rotations of the external momenta configuration. We consider only 
amplitudes of form 

\M\ 2 = \M\ 2 (s, Sl ) , (A3) 

where s = (px + Py) 2 an d s\ = [j>a — Vx) 2 - 

To obtain the differential decay rate, dT/ds, we must integrate out all other variables 
except for the Lorentz invariant s = (px +Py) 2 - To do this, we first integrate out py using 
the momentum conserving delta function. Changing to polar coordinates, we then have 

r = , L f \M^p^\P^h lM \ 2 5(m A -Ex-E Y - Ea) . (A4) 
2(27r) 5 m A y 8E X E Y E C K J 
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By Eq. (1A3j) . the choice for the axes with respect to which Qx and Qc are calculated is 
arbitrary. We choose the 3-axis to coincide with the direction of px so that 

dflc = d cos Qxc d(f>c , (A5) 

where 9xc is the angle between pc and px- Exploiting the simple relation 

EidEi = \pi\d\pi\ (A6) 

and performing the trivial dflx and d(f>c integrals, Eq. flA40 reduces to 

f 1 f \Px I \Pc I dE x dE c d cos Q XC ] K/n2sf w w P n / av ^ 

r = 7^3 / w \M\ 5{m A - E x - E Y - E c ) . A7 

{2-Ky i m A J 8E Y 

Note that the dcosOxc integral is non-trivial, since Ey is a function of cos^xc*- That is, 



e y = y \px + Pc\ 2 + m Y = sJm x +m 2 c + m Y + 2|p x ||pc| cos 9 X c ■ (A8) 

We are left with three integrations and one delta function. First, to conveniently encode 
the [—1, 1] domain of the cos 9xc integral, we add a theta function factor, 9(1 — cos Qxc 2 ) 
and extend the cos^xc integration limits to the entire real line. We now use the delta 
function to perform the integration over dcosOxc- From Eq. ( 1A8I) observe that 



dE Y _ \Px\\pc\ \Px\\p 



c\ 



d cos Qxc \/ m v + m c + m Y + 2|px||Pc| cosQ X c &y 



(A9) 



so that we may straightforwardly change variable from cos Qxc to Ey and perform the final 
delta function integral. We end up with 

T A ->x,y,c = / dE x dE c \M\ 2 Q{\ - cos 2 Q XC ) ■ (AlO) 

8(271)^711 A J 

Here cos Qxc is implicitly a function of Ex,c and the masses. 
We can now change variable from (E x , E c ) to (s, Si) using 



m A + m x ~ g i P m A + m c ~ s _^ it? dsds i 
Ztua 2rriA Am A 



then from Eq. (jAlOj) 



cis 32(27r) 3 m^ 



c?si|A^| 2 (s, si)0(l -cos 2 Q XC ) • (A12) 



In order to obtain the differential decay rate we need to perform just the si integration. 
This integral can be performed once an explicit expression for \Ai\ 2 is given. 

Before proceeding to consider the differential decay rates arising from different matrix 
elements, we still need to explicitly write down the limits of the S\ integration due to the 
9 function. In order to do so we must write cos Qxc as a function of (s, si) and then solve 
cos 2 Qxc < 1 for s i- Define s 2 = (px + Pc) 2 - Then 

s 2 = (px+Pc) 2 = m 2 x + m 2 c + 2E X E C - 2\p x \\p c \ cosQ XC • (A13) 
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Applying the identities (lAllj) . plus 

\Pi\ = 



EI -ml 



s + Si + s 2 = m 2 x + m 2 c + m\ + m\ , 



Eq. (1A13[) becomes 

2m A (s + S\ — my — m^) + (m^ + m\ — Si){m\ + — s) 



cos 9 



xc 



a/ (m 2 A + — s) 2 — Am 2 A m 2 c ^J (m 2 A + m 2 x — Si) 2 — Am\ 
The general solution for cos 2 8xc < 1 in terms of Si is S\_ < si < si with limits 



(A14) 



(A15) 



'i± 



m 2 A + m 2 x + my + m 2 c — s (my — m x )(m A — m 2 c ) 



2s 



± 



m a + m 



c 



m A m c 



2 2 i \ 2 

mf — m x + s 



— m Y s 



(A16) 



We now have our master formula for the differential decay rate 



dT 



ds 32(27r) 3 m 3 



A Js 



81, 



dsi\M\(s, si 



(A17) 



If the two leptons X and Y are massless, the integration limits reduce to 



m 2 A + m 2 c — s 



± 



2 i 2 \ 2 

+ — s 



— m^m 2 -, . 



(A18) 



Defining a new variable 77 = 1 — si/m 2 B , we immediately obtain the integration limits 
and the natural definition of £ as written in Eq. (fl2"j) of the main text. 

We now consider explicit examples and carry out the integral in (1A17j) for the three cases 
considered in the main text. 



1. (/> 3 interaction 

First consider the case in which all the particles involved are scalars with 3 vertices of 
the form 

gx4>A4>B4>x + 9y4>c4>b4>y • (A19) 

To take into account the decay width we use the Breit-Wigner approximation described 
in the main text ([7]). The matrix element is 

|A/fl 2 = gggj = 9 x9y (\nn\ 

' U ((p - Vx? - ml) 2 + m 2 B T 2 B - m|) 2 + m| Y\ ' 1 UJ 

so the master formula ( 1A17I) then becomes 

®- = gM f 1+ dSl 1 f A211 
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Using 



ds\- 



(si - m|) + m|r| m B r B 
we get the final result reported in the main text 

dT 



tan 



m B T 



9x9y 



tan 1 ( ^-r]. 
1 B 



ds 32(27r) 3 m^m s r B _ 
where % = 1 — si ± /m| as defined in the main text (cf. Eq. (|T2l ) . 



(A22) 



(A23) 



2. Intermediate Scalar 

We now move on to consider the decay in Fig. [Th. The corresponding Yukawa couplings 
are (Eq. ([32])) 



^ uk = M9lPl + 9RPRYxU ai ^ l B + M9l P L + 9lP R )^U ai *c 



B > 



(A24) 



which include flavor mixing. We denote the amplitude for the jth B mass eigenstate as 
Mf, and 



%Mf = u a (p x ) [g^P R +g^P L ]u(p A )u(p c ) [g Y L P L +g Y R P R y{p 



Y 



{p\ - m]) + irrijTj 

The total squared amplitude final involves both square and interference terms 
\M^ ot \ 2 = \MT\ 2 + \M S 2 C \ 2 + 2Re\MT(M s 2 c Y . 
We focus first on the interference term. After computing the traces, we have 



(A25) 



(A26) 



2Re 



Mf{Mff 



id? + (9rY islf + (g r R r MUaiU^u^ 



-X\2 



Y\2 



Y\2 



X 



[s 1 - m\)(s l - ml) + m 1 m 2 r 1 r 2 



s 1 — m\) 2 + m\T\ (si — ml) 2 + m|r 2 



7~ s ' 



and 



T sc = 64(p A ■ p x ){pc ■ Py) = 16(m 2 A - s 1 )(s 1 - m 2 c ) . 
Note that s\ = p\. The square terms |A4:- C | 2 follow from Eq. (1A25I) 



\M 



scl 2 



16 



{g x L ? + {g x R ? (gl) 2 + (g Y R ) 2 \u^\ 2 \u^\ 



[m\ - s 1 )(s 1 - m 



c) 



(si - m 2 ) 2 + m 2 r 2 



(A27) 
(A28) 
(A29) 

(A30) 



The differential decay rate can be obtained plugging Eqs. (1A27|) and (1A30P in the master 
formula (1A17I) . The final result is the sum of three integrations, involving the square terms 
and the interference one 



dT L 



ds 



{g x L ? + (g%) 2 } [{glf + (g Y R ) 2 ] 



(27r) 3 m3 



3=1,2 



12 



(A31) 
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Here, 



81, 



x j = \U aj \ 2 \U Pj \ 2 I ds v , 

lsi_ ( s l _ m j) + m j^j 



m 2 i \U aj \ 2 \U f3i \' 



+ — 



i + {-y z 



m 2 A + m 2 c 
2m^ 



4 + r 2 



m) 



tan 1 



rrij 

3 J 



1 lot 



|[/^| 2 |f/»|V 



ri + (-)v 


tan 1 


rrij 




. 3 . 











m. 



r]=rf + 



V=rf_ 



ri=7f + 



n=vL 
v=v J + 



r]=r]_ 



(A32) 



where m, z, y, and x are the oscillation parameters defined in Eqs. (J2"7]l and (128]) . r/^ = 

1 - s l± /m 2 , 



{m\ — m 2 )(m 2 — tuq) 



3J_\_ 3 
1 



(A33) 



and the approximation in Eq. ( 1A32I) drops terms whose prefactors are suppressed by Tj/rrij. 

The interference term is more involved. We report results to leading order in e and we 
directly write the coefficients in terms of the oscillation parameters as in the main text. We 
have 



X 12 = 2Re[U al m i U* 2 Up 2 ] / dsi 





(si - mf)(s 1 - ml) + m l m 2 T l T 2 




- m 2 ) 2 + m 2 Y 2 


(si — ml) 2 + m 2 2 Y 2 2 



RejU^U^U^U^m 2 
x 2 + (1 + yz) 2 



3=1,2 



r)=rf + 



V=V- 



>f I — 



77=77; 



(A34) 



where (at leading order in T/m = 2z/x) 

m 2 A m 2 (l — yz)(l — z 2 ) — m 2 A m 2 c (l + yz) + m 2 c m 2 (l — yz)(\ — z 2 ) 



A 3 



m H 



m 4 (l + {-) 3 zf[l + (-)M-3 + 3y + yz)] 



sc 2 m 4 



(A35) 



These coefficients and X 12 prefactors reduce to those in Eqs. fl38l) at leading order in z. 
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3. Intermediate Fermion 



Finally, suppose the intermediate particles are fermions. The Yukawa couplings are as in 
Eq. m 



4uk = {9lPl + 9i Pr) U ai *£^cf>\ + (g Y L P L + 9 Y r Pr) U m *i 
We denote the amplitude for the jth B mass eigenstate as Aij, and 



We ■ 



(A36) 



iM f = iU^U^*u a (p x ) [g x P R + g*P L ] 



B 



+ rrij — iTj/2 



bIPL + ^Pj^(py) • (A37) 



Now, since A and C are scalars, the total helicity of the \XY) final state must be zero. 
Since X and Y are massless, spin-1/2 particles, if they have the same (opposite) chirality then 
they must have opposite (same) sign components of momentum in the angular momentum 
direction. These are the only options. As a result, the terms of the differential decay rate 
may depend only on a linear combination of (g X g Y ) 2 , (9r9r) 2 , (.9r9l) 2 an d (9l9rY '' i an d 
the terms for the latter (former) two factors must have the same s dependence. Hence we 
may write the differential decay rate in the form (cf. Eqs. ( l4"3l) and ( 147|) ) 



ds 



M~9l) 2 + {9 X r~9 Y r?] 



dT a(3 



ds 



+ M~9 y r? + {g X R~9lf] 



dT a(3 



ds 



(A38) 



Here dT/ds\± are called the chiral decay rates for reasons explained in the main text, and 
each one of them has both square and interference contributions: 



dT a P 



ds 



dTl 



± 



ds 



dr 



aB 



ds 



+ 



± 



ds 



(A39) 



± 



For the sake of brevity, we introduce the following notation 

9 + = (9 X 9 Y r) 2 + {9 X R 9l)\ 9- = ^If + (g^lf 



(A40) 



and simply report here all the integrands in terms of the variables of integration (s, si). The 
interference term is 



M { 



32g + Re{U al m i U* a2 U l 



12 



[(si - m 2 ) 2 + m?rf] [(si - m|) 2 + m 2 2 T 2 2 ] 



x < mxm 2 s 



(si - m\)(si - ml) + m 1 m 2 T 1 T 2 



1 

+ 2 



si - m\)m{Ti - (s\ - mj)m 2 T 2 (miT 2 - m^^s 
32g_Re{U al U*p 1 U* 2 Up 2 ] 



(A41) 



12 



[(si - m 2 ) 2 + mfrj] [(si - ml) 2 + m 2 T 2 ] 

x [{m\ - si)(s 1 - 771^) ~~ ~~ m l)( s i - m i) + mim 2 r 1 r 2 ] . (A42) 
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In contrast, the square terms are 

fl9 32 g + \U^\ 2 \U^\ 2 m 2 s 
M i l = u 2 2 ' 2r i > (A43) 

I w£, a = 32 ff-|^| 2 |^| 2 [(mj - giKgi - ml) - ssj 

1 iU [( ai -mJ)2 + mjrj] ' 1 ' 

In order to obtain the final expression for dT a ^ /dsL we just need to perform the integra- 
tions over si as in Eq. ( ]A17j) . At leading order in e, the results are displayed in Eqs. ( I4"9j) 
and (I5UI) . but the full coefficients of the interference terms (including terms to all z orders) 
are 

A j - — 

2 m z 
A j = A 1 - — 

• x sm 2 

B - + (-) 5 2^f " ( A45 ) 

with ^ c and ^4| c defined in Eq. (1A35I1 . Note that in the vectorial coupling case, the linear 
s dependence of the log coefficient terms cancel up to terms of order 



xz x 2 



X 2 + 1 X 2 + 1 

which are negligible, as expected. 
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